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, ABSTRACT. Just like decent classical difference-difference systems 

define symplectic maps on suitable phase spaces, their counterparts with 
properly ordered noncommutative entries come as Heisenberg equations 
, of motion for corresponding quantum discrete-discrete models. We ob- 

serve how this idea applies to a difference-difference counterpart of the 
Liouville equation. We produce explicit forms of of its evolution oper- 
CN ■ ator for the two natural space-time coordinate systems. We discover 

Jjf" ' that discrete-discrete models inherit crucial features of their continuous- 

^r) ' time parents like locality and integrability while the new-found algebraic 

CD ■ transparency promises a useful progress in some branches of Quantum 

Inverse Scattering Method. 
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^ ■ 1 Introduction 

<D ' 

■ In this paper we intend to elucidate the algebraic part of our approach to the quan- 

tum integrable models in 1+1 dimensional discrete space-time, developed during 
last five years. We shall not give a complete survey of our publications (Faddeev 



X 



5_j ■ and Volkov 1993,1994; Faddeev 1994; Volkov 1997a,b) because it would take too 



C3 



much space. We believe, that the algebraic side is most instructive and original; 
more analytic side will be mentioned only briefly with references to (Faddeev 1994) 

Discrete space-time models (DSTM) in soliton theory have acquired a prominent 
role from the very advent of this part of mathematical physics. The first examples 
of such models were proposed by Hirota (1977) as discrete analogues of the major 
continuous soliton models. Subsequent development was carried on mostly by Dutch 
group, see (Nijhoff and Capel 1995) and references therein. The recent resurgence 
of interest towards DSTM is connected with several new ideas: 

1. The nonlinear equations for the family of transfer matrices Is (A) in the 
framework of the Thermodynamic Bethe Ansatz can be considered as DSTM with 
spin S and rapidity A being discrete variables (Kliimper and Pearce 1992; Kuniba 
et al. 1994; Krichever et al. 1996). Of course, the rapidity assumes the continuous 
values, but only discrete shift A — > qX enters the equations. 
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2. A. Bobenko and U. Pinkal with collaborators develop the discrete analogue 
of classical continuous 2-dimensional differential geometry, see (Bobenko and Pinkal 
199?) and references therein. 

3. Quantum version of DSTM revealed a new type of symmetry, giving the 
discrete analogue of current algebra and Virasoro algebra (Faddeev and Volkov 1993; 
Volkov 1997c). Moreover, quantum DSTM seems to be rather universal, giving both 
massless (conformal) and massive models in continuous limit. 

For the evident methodological reason we shall illustrate our approach on a 
concrete example. For that we have chosen the most prominent model of physics 
and geometry — the Liouville model. More involved Sine-Gordon model will be 
touched upon only briefly. Incidentally, the latter was already a subject of our 
earlier publications. We shall not discuss the usual paraphernalia of integrable 
models such as zero-curvature representation, Lax equation and Bethe ansatz. We 
shall simply present the main dynamical object — the evolution operator, realizing 
the elementary time-shift. Its natural place inside the Algebraic Bethe Ansatz is 
discussed in recent lectures of one of authors (Faddeev 1996). We believe, that our 
explicit formulas are interesting enough as they stand so we want to present them 
in its clearest form, independent of original derivation. 

We begin with the reminder of the classical Liouville model and its Hamiltonian 
interpretation. This will play the role of the starting point for the subsequent 
deformations: discretization of space on which hamiltonian data are given, and 
quantization. As a result we shall get a suitable algebra of observables. Finally the 
time evolution will be defined in terms of a certain automorphism of this algebra. 
The discrete time equations of motion produced by this automorphism will be shown 
to be a natural analogue of the corresponding classical equations. The integrability 
of the model will be confirmed by presenting an explicit set of conservation laws. 

2 Classical differential equation 

As the goals declared in the Introduction suggest, this time we shall consider the 
evergreen Liouville's Equation (LE) 



a Hamiltonian 1+1-dimensional field theory with x denoting the spatial coordinate 
and t serving as time. The Cauchy data 



8P<p d 2 ip 



dt 2 dx 2 



<p(x,t)\ 



t=0 



tp{x) 




t=o 



w{x) 



can be equipped with the canonical Poisson bracket 



{w(x),ip(y)} = 5{x-y) 



{w(x),w(y)} = {ip(x),ip(y)} = 
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so that the evolution goes the Hamiltonian way in its most familiar 
if = {H, (p} = w w = {H, tu} = if" + e -2 ^, 

the Hamiltonian being 

H = | J dx (w 2 + (ip') 2 + e-^). 

The periodic boundary conditions 

ip(x + 2tt) = <f(x) w(x + 27r) = w(x) 

pose no problem provided the Poisson bracket employs a 27r-periodic delta-function 
rather than the ordinary one. 

Of the equation's specific features the Liouville's formula 

(f(0-9(r)) 2 

x = ^ — t t = £ + r 

making a solution out of two arbitrary functions, is the ultimate. It is there, ac- 
cording to (Gervais and Neveu 1982; Faddeev and Takhtajan 1985), where the real 
Hamiltonian theory of LE begins. We shall not reach that high in this paper. 

3 Classical difference equation 

The best lattice approximation to LE 

e ip(x,t-A) e <p(x,t+A) _ eV >(x-A,t) e tp(x+A,t) = A 2 

is due to R. Hirota (1987) like virtually every decent difference-difference equation. 
In order to make its transformation into LE under limit A — > more obvious one 
may recompose it like this: 

sinh | (^p{x, t - A) + <p{x, t + A) - <p(x - A, t) - fix + A, t)j 

= /\2 e -^Mx,t-A)+<p(x,t+A)+<p(x-A,t)+<p(x+A,t)) _ 

Now as the mission of the lattice spacing A is over, it is only natural to have 
everything suitably rescaled 

(x,t) — ► (Ax, At) 
e v — > Ae^ 
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or just set A=l. Either way, the Difference Liouville Equation (DLE) takes its final 
form 

gVj.fc+igVj.fc-i _ gVj+i.fcgVj-i,* = x 
j + k = 1 (mod 2) 

where the change for subscripts manifests that the 'space-time' is now a Z 2 lattice 
while the second line specifies which half of that lattice the equation will actually 
occupy. This half itself makes a square lattice turned by fourty five degrees with 
respect to the original one and twice less dense. The values of 92 on a 'saw' formed 
by vertices with k equal either or 1 

fj,o = <Pj for even j 

(p jt i = ipj for odd j 

make a reasonable Cauchy data, that is they are just sufficient to have the whole 
system resolved step by step. This has everything to do with the second-order nature 
of the original continuum equation whose Cauchy data combine the present and a 
little bit of the future represented by <p(x) and w(x) respectively. 

Quite expectedly, there exists a Poisson bracket preserved under evolution along 
/c-direction governed by DLE. However, it turns out more complicated than one 
might have wished a lattice deformation of the canonical one would be: 

{<Pi,<Pj} = 

with 

{0 if i - j = 1 (mod 2) 

(—1) 2 ( i+:,+1 )sign(i — j) otherwise 

Such is the price for the ultimate simplicity of the equation. This would be too 
much if we had lost the option of periodic boundary condition 

Pj+L = <Pj. 

Fortunately, we had not. If the period is chosen properly 

L = 2M M = 1 (mod 2) 

the bracket remains intact provided the above description of ? applies when — < 
L and extends periodically 

elsewhere. Those still insisting on an easier bracket can change the variables 

<A? = hiw+i + ¥>j-i) 
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and have it 

{^■-i.^-} = 1 Y L 

at the expence of a busier equation 

^.4>j,k+i e ^<l>i,k-i = (1 + e 2 ^' +1 '*)(l + e 2< ^' _1 ' fc ) 
j + A; = (mod 2). 

Either way, the prospect of dealing with discrete Poisson maps is hardly encouraging. 
That is why we choose to leave the classical equations alone and go quantum. Before 
we do, let us round out the classical part with a beautiful, if irrelevant for our current 
agenda, discrete Liouville formula: 

e -20 jife = e -V j+ i,k-Vj-i,k = (fm+1 ~ /m)(ffn+l ~ 9n) 

(fm+1 — 9n)(fm ~ 9n+l) 

j = m — n k = m + n + 1. 

4 Algebra of observables 

One dilemma about quantization is whether to develop it in terms of the bare <p's 
or stick to the variables actually entering the equation, that is the exponents 

vj = e Vi {vi, vj} = q(i,j)viVj 

The respective Heisenberg- and Weyl-style quantum algebras are 

[<Pi,<Pj] = 

<Pj+L = ¥>, 

with the usual lot in r.h.s. comprising the imaginary unit i, the Plank constant h 
and the coupling constant 7; and 

V j+L = Vj 

with all packed in a single quantisation constant 

q = e lh \ 
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Weyl-style algebra may be viewed as a subalgebra of the Heisenberg-one 

Xtj = e*i 

but not the other way round. Roughly speaking, the latter also accomodates non- 
integer powers 

« = e a ^ 

not allowed in the former. 

Another dilemma is whether to place q on the unit circle or not. The first option 
is obviously involution-friendly. It allows for both unitary 

and self adjoint 

»i = »i 

pictures, the former offering the luxury of dealing with bounded operators if at the 
expence of complications in representation theory due to the arithmetics of q while 
the latter actually being the one relevant for the true Liouville model. On the other 
hand, q inside (or outside) the circle is favoured in q-algebra but whether it is good 
for something else remains to be seen. 

We choose not to take sides before time and conclude the Section on a more 
practical note. Let us introduce, for future use, quantum counterparts of the e 2 ^- 
variables 

tOj = Uj + it>j_i = 0j_i0j + i 
and compile a list of emerging commutation relations 

tOj-itOj = q 2 ^ 1 ^ 3 K>jtt>j-i 

ttitVj = tOjVi if % ^ j (mod L) 

tOiKij = KijKii if \i — j\ ^ 1 (mod L). 

5 Evolution operator 

Given an invertible operator 0, one can make the algebra of observables 'evolve' 

• • • i — ► QiQ.' 1 1 — > 3 1 — ► 1 — > n _2 3£i 2 1 — > ■■■ 
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hoping that the evolution of generators 

iife+2 = Q _1 t)j, fe £3 j + k = (mod 2) 

l>2a,0 = 02a f2a-l,l = t>2a-l 

manages to solve some nice and local equations, for instance, 

t> j,k-i ~ 9 _1 t>j-i,fcf = !• 
As a matter of fact, this is exactly what happens if 

M M 

£2 = IJe(m 2 o-i) 5 ri 6 ^) 

a=l a=l 

provided 5 is the 'flip' operator 

r V = "7 1 

while the function e solves the following functional equation: 

f(gg) _ 1 
e(q~ 1 z) 1 + z 

Indeed, let us plug the definition of the tJj^'s into the hypothetical equation 

Q- b - 1 t> 2a Q b+1 Q-%aQ b ~ q-'Q-^a-l&Q-^a+lQ' = 1 
£T b 9 2a -lCl b £r b+1 X>2a-l£l b - 1 ~ q- 1 Q-%a-2& b n- b V2aQ b = 1 

and dispose of as many 0's as possible: 

O2a0O 2a - q~ 1 Q-V2a-lt>2a+l = H 
02a-in0 2a -l - g _1 t>2a-2t>2a£2 = £2 

Then all the e(tt>) factors but one go the same way which results in 
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Once 5 is gone too, we are left with 



°2a e ("'2a)f 2a ~ Q ^2a)^2a-1^2a+l = e («'2a) 
*>2a~AK>2a~l)K2a-l ~ 1~ lv 2a-2»2at(K>2a-l) = ^2a-l) 

which is nothing but the above functional equation mated with the commutation 
relations which closed the last Section: 

»2a e (K>2a>2a = <Q^2a) ^>2a-A^2a-l>2a-l = <Q^2a-l) 

e(q- 2 t0j) - q-^itorftOj = e(q- 2 t0j) - q^tv^Wj) = e(tOj) . 

So, since everything eventually reduces to that functional equation, let us see if it 
can be solved. 

6 q- exponent 

Indeed, the equation in question 

<qz) _ 1 
is readily fulfilled by those ubiquitous q-exponents 

e(*0 = (-Q z ;q 2 )oo - good for |g| < 1 

£ ( Z ) = 7 -1 1 _9N - g°° d fOT M > 1 

(-q l z;q 

where 

oo 

There is no solution entire in z if \q\ = 1. For those opted for the Weyl-style algebra 
of observables (see Section 3) this is the end to the story, not a happy one in the 
latter case where the equations of motion survive but the evolution automorphism 
behind them turns outer. 

The Heisenberg Way has more solutions to its disposal: 

eHeisenbcrg(^) = eweyl(z) x any function^) 
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with 9 being the first solution 

to the equation 
with a clear purpose: 



7T 

ftrf 



q 26 = l 



Among them we find the one capable of surviving under \q\ — > 1 limit (Faddeev 
1994): 

(-^;? 2 )oo 



e g (z) = e*(z e ) = < 



(-?-»V; 
(-q $2 z 9 ;q 2e2 ) 



\q\ < 1 



I (-^ 2 ) c 

It is plain to see what makes it so special. Duality is the word. We get two functional 
equations 



e q (qz) 1 
e g (g -1 2) 1 + z 



satisfied at once. Consequently, the updated evolution operator 

M MM M 

£ = II ^(tXJ 2a _!) 5 n e ?( W 2a) = II V 2 ^-l) * II V 2 ^a) 



o=l 



a=l 



a=l 



o=l 



not only inherits the right evolution of the pure d's 

t>j,fc+lt>j,fc_l - ^"^j-l^tij + l,*; = 1 

but also produces an equally right evolution 

"ifc+i^ife-i ~ 1 6 °j-i,kQ e j+i,k = 1 

of their dual twins s . Since the Heisenberg-setup turned out to be just a pair of 
decoupled Weyl-ones, we will stick to the latter for the rest of the paper. 
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7 A different angle 



Although the ingredients used in the formula of the evolution operator are all more 
or less familiar, they are put together in a bizarre way. A traditional R-matrix 
philosophy would offer a different approach which we now start presenting. The 
equation remains the same 

v m+l,n+l v ra,n — Q 1 ^m,n+l^m+l,n = 1 

but the change for a kind of 'light-cone' variables 

j = m — n k = m + n 

signals that the n-direction is now considered temporal. So, we are going to find 
out what algebra of observables and what evolution operator produce the evolution 

V m ,n+1 = Q lv m,nQ 
Vm,0 — V m 

solving that 'light-cone' equation. 

8 Algebra of observables ii 

Here is the complete list of relations defining the new algebra of observables: 

y £ y m = y m y £ if m - £ = 0, 2, . . . , M - 1 

V£V m = q\i m vi if m = 1, 3, . . . , M 
v m C = qCv m 

v m +M = <?^Cv m M = 1 (mod 2). 

This is exactly what it takes to achieve the required relationship 

v m w m = q 2 w m y m 

V£W m = w m V£ if £ ^ m (mod M) 



10 



between the v's and their 'derivatives' 

v m+l 

w m = 

Vm-l 

which themselves form the much advertized by the authors 'lattice current algebra' 

2 

w m -iw m = q w m w m _i 

\N(\N m = w m W£ if \m — l\ 7^ 1 (mod M) 
W m +Af = w m . 



We already met similar relations, it was the end of Section 3. That time they did 
not contradict the periodicity of the D's. Now they do: it is impossible to have C = 1 
and good commutation relations at the same time. We shall soon see why. 

9 Evolution operator ii 

Let us see what the operator 

Qnaive = SFe(wM) • • • e(w 2 )e(wi) 

can do. Of course, the function e and the flip operator F are the same as before 

e(gz) = 1 
e(g~ 1 z) 1 + z 

F _1 v F = v _1 

while S is the shift operator 

S 1 v m S = v m _i 
We plug the 'naive' evolution into the hypothetical equation: 

Q" n_1 v m+ iQ n+1 Q" n v m Q n ■ - o- 1 Q" n_1 v m Q n+1 Q" n v m +iQ n ■ =1 

naive v ™+± ^ naivc^* naive mU < naive V ^ naive '"^'naivc^'naivc m + lu < naive 

and dispose of as many Q na ive' s as possible: 

v m+lQnaivc v m Q v m.Qnaivc v m+l = Qnaive- 

Then all the e(w) factors but one go the same way which results in 
v m+ iSFe(w m )v m - g~ 1 v m SFe(w m )v m+ i = SFe(w m ). 
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Once S and F are gone too, we are left with 

v m le ( W m) V m ~ ?~ lv m-l e ( W m) v m+l = e( w m) 

which is nothing but our functional equation mated with the commutation relations 
which closed the last Section: 

v m 1 e(w m )v m = e(<T 2 w m ) 

e(<T 2 w m ) - ^ _1 e(w m )w m = e(w m ). 

This proves that the 'naive' evolution satisfies the required equations of motion ... 
as long as m is neither 1 nor M. We could not reasonably expect it to do any better 
because Q na ive obviously had no respect to the translational symmetry of the algebra 
of observables. In order to have this eventually repaired, let us first figure out how 
that sad dependence on the starting point can be cured in a simpler situation, say, 
for a monomial . . . w^w^ 1 . Pulling the factors from the very right to the very 
left one by one we get a clear picture of how the ordered monomials with matching 
powers but different starting points turn into each other: 

w^...wfwf 

= q 2 PMPl q -2 P lP2 W P1 W PM _ W |3 W P2 
= q 2 PMPl q -2 P 2P 3 W P2 W P1 ^4^4 = 

Now we know. The expression q- 2 P™-iPm w^S^w^Z^ ■ • • w ^+i does not depend 
on m provided pe+M = Pe- We award it with a self-explanatory notation 

O 

11 wf = q- 2 PMPi w Pm _ _ w P2 w Pi 
and extend this definition linearly to the corresponding polynomials, in particular, 

n e ( w ^) = Yj c pm • • • °p2 c pi n wPe 

PM, — ,P2,P1 

r ■■■r r a~ 2 P MPl \n Pm w P2 w Pl 
~~ P M c P2 C Pi1 W M ■•• w 2 w l 

PM,— ,P2,P1 

= £ C pM C pi g- 2 ^ 1 W^(' e (w M _ 1 )... e (w3) e (w 2 ))wf 
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- E 

Pm — l iPn 



l C PmQ 



-2p m _lPn 



w 



Pm-1 

m— 1 



e(w. 



m-2; 



••e(w m+2 )e(v m+ i) 



with the coefficients c coming from 

p 

The same treatment applies as well to the 'selfdual' e's of Section 5: 

P,r 

PM,Pi,rM,ri 

xw m /w m M ( e <?( w M-i) • • • e 9 (w 3 )e 9 (w 2 )^wf wf 
So, we achieve the vital translational invariance of those products 



sn e ( w ^)=ii e ( w ^ s 

sacrificing none of their 'orderness'. The time has come to plug the repaired evolu- 
tion operator 



Q = SFj]e(w,) 

into the hypothetical equation ... see the beginning of this Section. 

Now as we finally established that the operator Q is indeed responsible for the 
quantized and fully discretized Liouville equation 

V m +l,n+l v m,ri — q 1 V m ,n+l V m +l,n = 1> 

we must admit that so far the commitment to this particular equation was only a 
matter of personal taste. What would change if the function involved 

O 

Q = SFj]/(w,) 

was different, for instance, 

/(*)- €{Z) ' 



e(q 2X z) 
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Nothing except the r.h.s. of the functional equation 

f(qz) 1 + q 2X z 

and the form of the eventual equations of motion 

V m +l,n+l v m,n ~~ 1 lv m,n+l v m+l,n = 1 — q X+1 ^m+l,n+l^m,n+l^m+l,n^m,n- 

By the way, this is another Hirota's equation, the discrete sine-Gordon one. We 
shall see it again in Section 11. 

10 Classical continuum limit 

The equation of Section 6 certainly turns into the 'light-cone' Liouville equation 

didr 

the matching Cauchy problem being 

y>(£,T)| r=0 =y(0- 

The algebra of observables from Section 7 transforms into no less familiar Poisson 
bracket reading 

{tp(0,tp(v)} = isign(£-r/). 

The evolution operator of Section 8 bears some resemblance to the corresponding 
Hamiltonian 

What is wrong? We seem to inherit also the quasiperiodic boundary condition 

which obviously contradicts to the equation. The periodic condition 

^ + 7 r) = Y>(0 

could do but that in turn would contradict the Poisson bracket. A more careful 
examination reveals that the 'constant' in boundary conditions is not a constant of 
motion: 

Q _1 CQ = C _1 . 

Of course, the lattice equations of motion themselves have no problem with that. 
However, their solutions are not smooth enough to survive a straightforward con- 
tinuum limit. Anyway, this peculiarity is not too relevant to what we are after in 
this paper. 
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11 Conservation laws 



Looking at the two evolution operators we now possess 

MM 
H = [J e (">2a-l) d J] e ("'2a) Q = SF f[ e( W< ) 

a=l a=l 

do we see something in the latter that was not there in the former? We do, the latter 
looks almost like a good old ordered product of 'fundamental R-matrices' (Tarasov 
et al. 1983). According to (Volkov 1997a), the shift-n-flipless part of Q 

O 

Q = JJe(w^) 

belongs 
in a family 

consolidated by the Artin- Yang-Baxter's Equation 

Rm+ltVjRmWRm+lM = R m MRm+lWR m (^f)- 

The choice of notation 

Rm(A) 



Rm(A) = e(A|w m ) Rm(A,M) = 



RmW 



is meant to emphasize the R-matrix connection. Let us recall how that AYBE could 
be verified. From (Faddeev and Volkov 1993) comes the multiplication rule 

e(A|b)e(A|a) = e(A|a + b + gba) 

applying whenever a and b satisfy the Weyl's algebra 

ab = (/ 2 ba. 

The two w's next to each other certainly do, therefore 

R m+ i(A)R m (A) = e(A|w m + w m+ i + gw m+:L w m ), 

therefore 

Rm+l (A) R m (A) R m +l (m) Rm (m) 
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= e(A|w m + w m+ i + gw m+ iw m )e( / t/|w m + w m+ i + gw m+ iw m ) 
= e(ji\w m + w m+1 + gw m+1 w m )e(A|w m + w m+ i + qvj m +i^ m ) 

Rm+1 O) Rm (m) Rm+1 (A) R m (A) . 



This is it. 

AYBE and ordered products are known to make a natural match 



RmW^m-lW • • • Rl( A ) Rm+lWRmM • • • R2M 



Rm(A)R m +l(/i) ) ( R m— 1 

(A)R m ( M ) ... Ri(a)R 2 (m) 



Rm+l(A) 

_ Rm+1 (A) 

= Rm(M)Rm+l (A))j^^R m—1 

(A)R m ( M )J ... ( Ri(a)R 2 (m) 

= ^R,„-/'-R,„ . I A ^ ^R„, p//-R ( „.A^ I . . j RjfAjR^f//) 

= ^RmMRm+lwj ^Rm-lWRm(A)j ■ ■ ■ ^ j ( R| (A)Ro(/') 



(A) R m—1 (M)Rm(A) ) ... ( Rl(M)R 2 (A) jRl(A, M ) 

Rl(A) 



RmWRm-lW • • • Rl(A> Rm+l(A)R m (A) • . . R 2 (A) 



RiM ! 



so, one hopes the O-ed product to step beyond m = M — 2 and deliver 

fi(A)fi(/x) = fi(/i)0(A). 

Let us see. 

fi(A)fi(/x) 

E r (A) „(A) -2p M pi w PM p(A) r( A ) w P! 
C Pm C Pi ^ W M R M-1 • • • K 2 w l 
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E r 0) r 0) n -2r M ri r M p (m) R^ w r i 
c r M c n y W M n M-l ■ ■ ■ n 2 W l 

E r ( x ) M r ( x ) M n -2(PMPi+rMri+r M Pi) 



xw PM R (A) w rM R (A) R (/l) R (A) R (At) w Pl R (M) w ri 
XW M R M-1 W M R M-2 r V-l • • • R 2 R 3 W 1 R 2 w l 



- so far, we only recalled the definitions and did some reshuffling not involving any 
nontrivial commutation relations, (A) and moved to the superscript level in order 
to save some paper - 



- r( A ) A^r M n - 2 (PMPl+rMri+r M pi) 
PM U r M ^pi v ri y 



d(p) d( a ) 



xw m K M-i w Af rTvrrTvr k m-2 k m-i ••• K 2 K 3 w l K 2 w l 



pO) rM 



the unit operator -^--^y 3 - = 1 has been inserted 

R M-1 R M-1 



- V r( A ) f (M) r ( A ) f (M)„-2(pMPl+»*Am+rMPl) 

R M r(A) 
v»a,P m r( A ) ,./m M-l D (p) R (A) d(m) d (A) k 2 ,J1dM,,Ji 
XW M K M-1 W M -(A)~ K M-2 K M-1 • • • K 2 K 3 T^Tj" W l K 2 W l 

°M-1 R 2 



- AYBE did its habitual job 



E r W M r W M n -2(pMPi+rMri+r M pi) 
^pm^tm^pi °n y 



d(m) p (A) 

XW M K M-1 W M (A) • • • r Sji+l r Sre+2 r Sn VlVrm ■■■ p (/i) 1 2 W l 
R M-1 R 2 



- nothing happened, we just refocused the attention to the middle portion of the 
product - 
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E r W>)>) M r W r W r WM n -2(PMPi+r M ri+r M Pi) 
Pm I'm V+iW+iV-pi °ri « 

R (f) R (A) 
via/P m rM ... r M M-l r m+ i R (A) Jm^m+lpW ...Pm ^2 ,.,PipO),.,n 

XW M K M-1 W M T(A)- • • • W m+1 K m+2 W m W m+1 K m-l W m ■ ■ ■ r^ W l K 2 W l 

R M-1 R 2 



- we disassembled some of the R's - 



r f (^) f «>) r (A) >) r (A)„(A) („) 



p(m) 

W M K M-1 W M "(AT • • • W m+1 K m+2 W m+l 
R M-1 

(A) 



R 2 



- the two in the middle traded places 



- \^ JaOJA) JA) ( M ) (A) (ji) Jji) (A) -2(r m r m+ i+p m p m+ i+p m r m+ i-rip M ) 
/ R (A) n 

x(w-R^ 1 w^...^ y wfR^w-j 
/ R^ 

X ( W M K M-1 W M "(A) - • • • W m+1 R m+2 W m+l 
V R M-1 

- the two halves in big brackets passed through each other - 

-= » rW rW » M M A X ) 

Z-> I'm V Pi ^PM^l ^M^m+l Pm+I y 

R (A) r (m) 
xw f™RW w Pm ^2 Pi p (m) w Pm w ri R W r M M-l w r ™+i p( A ) W P™+1 
XW m R m-l W m ■•■ (^) w l K 2 W M W l R M-1 W M (A) 1 ' ' W m+1 R m+2 W m+l 

R 2 R M-1 



- the two in the middle traded places 
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E 



c M c ( a ) c (m) c (A) -2(r m r m+ i+p m p m+1 +p m r m+1 ) 



D (A) dM 



Xw m K m-l w m ■ ■ -T(^)" K 1 K 2 K M K l K M-1 K M "(A) - - ■ - W m+1 K m+2 W m+l 
R 2 R Af-l 

- we assembled some R's, now it only remains to apply AYBE three more times 



— \^M r WM r W „-2(r m r m+ i+p m p m+ i+p m r m+ i) 

xw rm R (M) w Pm rMpWrWrWrM p>) W rm+1 R (A) w Pm+1 
xw m t\ m _jW m . . . t\ x r\ 2 r\ M r\ x r^^r^ • • • ™ m+1 ^ m+2 w m+l 

xVr (A) r (A) „- 2 P-P™+i w Pra R (A) R (A) w Pm+1 
= ^Ou)fi(A). 

Done, at least for M > 5. In fact, even M = 3 is possible but this would take 
three more pages to verify. Anyway, a more civilized edition of the above proof is 
presented in (Volkov 1997b). 

The commutativity of the f2's may be a good news but there is a bad one too. 
The flip operator F does not commute with the O's. Which means there is another 
family 

U(\) = F _1 fi(A)F = FSlMF- 1 
not coinciding with the original one. Of course, the ?J's commute with each other 

U(A)U(aO = U(aOU(a) 

but it is not immediately clear whether 

Q(\)l3(ii) = U(/i)0(A) 

should also be true. Fortunately, there is some hidden agenda making it happen. 
Technically-wise, the proof is the same as that above except the AYBE in use is 
somewhat different: 



( W m+l e ( A -^l W m+l)) e(A|w m )e(/i|w. 



m+1 J 
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= e(A-/i|w m+1 )e(A|g 2At w m )w^ +1 e(^|w m 1 +1 ) 



e((i\q- 2 »\N m+1 ) 



^^(H^w Je( A|^w m+ 0^™ 
= w^e(^|w m 1 )e(A|g- 2 ^w m+1 )e(A-^[w m ) 

= e(M|w~ 1 )e(A|w m+1 ) ( e(A-^|w m )w^ ) . 



Once the flip-n-shift join 

Q(A) = SFft(A) = ?J(A)SF 
one realizes that full commutativity is not there 

Q(A)QM + QwQ(A), 

only 'squares' can do: 



Q(k)Q(a)J |^Q(^)Q»J = FQ( K )FQ{\)FQ(fi)FQ(u) 
= F 4 c3( K )$7(A)t3( M )^H = F 4 0(^)^Hc3(k)0(a) 
= Fn(/i)FO(i/)FO( K )FO(A) = ( Q(p)Q(u)] (q(k)Q(X) 

In particular, 



CT 2 Q 2 (A)Q 2 = Q 2 (A). 

On these grounds, let us call the Q(A)'s 'conservation laws' even though what actually 
happens is that only their squares are only recovered on every other step in time. 
This peculiarity has everything to do with that discussed in Section 9. Apparently, 
not one but two time steps should make a 'physical' time unit. 

12 Conclusion 

We developed here the scheme allowing to describe some quantum dynamical sys- 
tems in discrete 1+1-dimensional space-time. The discretized Liouville model was 
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taken as the main example and treated both for laboratory coordinates and light- 
like ones. All considerations were purely algebraic, no representation and/or Hilbert 
space was used. We confined ourselves to pure Heisenberg picture of quantum the- 
ory. 

The main outcome is the construction of evolution operator realizing the au- 
tomorphism of the algebra of observables leading to the Heisenberg equations of 
motion representing lattice and quantum deformation of the corresponding classical 
equations. In this construction the famous q-exponent (q-dilogarithm) played the 
most prominent part. 

We discussed also the integrability of the model presenting the set of conservation 
laws. Their construction and the verification of commutativity was based on a new 
solution of Artin- Yang-Baxter relation, being a close relative of the q-exponent. 

We hope that the scheme of this paper is general enough and allows to include 
many related models of quantum theory. Our papers mentioned in Introduction 
give some illustration of this. 
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